CONVERSION OF MAGNETIC FIELD ENERGY INTO 

KINETIC ENERGY IN THE SOLAR WIND Y. C. Whang 


ABSTRACT The outflow of the solar magnetic field energy (the radial component of the Poynting 
vector) per steradian is inversely proportional to the solar wind velocity. It is a decreasing 
function of the heliocentric distance. When the magnetic field effect is included in the 
one-fluid model of the solar wind, the transformation of magnetic field energy into 
kinetic energy during the expansion process increases the solar wind velocity at 1 AU by 
17 percent. The predicted solar wind conditions at 1 AU are w = 302 km/sec, 
M “^protons/cm^, T- 1.5X10^°K, 1.4X10“^ erg/cm^/sec,^ = 7.37 0 = 129.5°, and 

P = 1.58. They agree very well with the observed quiet solar wind. 


mTRODUCTION 

Since the original application of hydrodynamic equa- 
tions to the solar wind by Parker [1958] , a number of 
theoretical models have been developed to study the 
solar wind, notably the one-fluid models studied by 
Noble and Scarf [1963] and by Whang and Chang 
[1965 ] , and the two-fluid model by S tarrock and Hartle 
[1966] and by Hartle and Sturrock [1968]. The 
one-fluid model assumed that all types of particles 
(electrons, proton, helium, etc.) have the same 
temperature, while the two-fluid model took into 
consideration the different temperatures for protons and 
electrons. The inviscid model of Whang and Chang 
predicted a solar wind velocity of 260 km/sec at 1 AU. 
The solution of Hartle and Sturrock predicted a velocity 
of 250 km/sec at 1 AU. These velocites are about 20 
percent below the observed quiet solar wind condition 
[Hundhausen, 1970] . 

The interplanetary magnetic field has a strong coupling 
with the dynamics of the solar wind. The magnetic field 
organizes the collisionless plasma into a continuum in 
the expansion process of the solar wind. It is well known 
that when the solar wind flows across a shock wave, the 
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magnetic field energy flux and the convective thermal 
energy increase during the compression process at the 
expense of the kinetic energy. In a reverse process, 
during the magnetohydro dynamic expansion of the solar 
wind, the magnetic field energy as well as the thermal 
energy is continuously converted into kinetic energy. 
Therefore, we propose that in a magnetohydro dynamic 
model of the solar wind, which includes an additional 
energy source, the predicted solar wind velocity can be 
increased closer to the observed quiet condition. 

Various mechanisms for the additional energy source 
needed to increase the solar wind velocity have been 
proposed by many authors. Barnes [1968, 1969] 
proposed that a significant amount of energy is carried 
by the hydromagnetic waves. The dissipation of 
hydromagnetic waves adds its energy into the solar wind 
flow. Jokipii and Davis [1969] proposed that the 
observed long-lived velocity streams are the source of 
long-wavelength turbulence, which in turn serves as an 
efficient and variable heat source for the solar plasma. In 
this paper, we will show that the magnetic field energy 
flow per steradian is a decreasing function of the 
heliocentric distance r. The total decrease of field energy 
flow from r = 2rQ to r - 1 AU is about 25 percent of the 
convective kinetic energy at 1 AU. Therefore, during the 
expansion process of the solar wind, the conversion of 
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field energy into kinetic energy increases the solar wind 
velocity by 17 percent. 

The solution of a steady-state magnetohydrodynamic 
one-fluid model of the solar wind is carried out in this 
paper. The model assumes that the solar wind velocity is 
in radial directions away from the center of the sun, and 
its temperature is isotropic. The solar wind conditions 
near the equatorial plane predicted by the present 
solution clearly demonstrate the process of energy 
conversion in the solar wind on a quantitative basis. 

DECREASE OF THE MAGNETIC FIELD 
ENERGY FLOW 

Spiral Pattern of the Magnetic Field 

The steady-state solar magnetic field is governed by 
Maxwell’s equations 


V-B = 0 

(1) 

V X(uX B) - 0 

(2) 


Spherical coordinates (/*, i//, co) are used in this paper; ij/ 
denotes the colatitude measured from the north pole of 
the sun. If B = e^Bf^ry\p) + then we can 

write equations (1) and (2) as 

i - 0 

These two equations can be integrated to give 

Br<=^r^ (3) 

B^^(ru)-^ (4) 

which represent the spiral pattern of the solar magnetic 
field. 

The Poynting Vector 

The equation of energy conservation can be integrated 
once: 

/ miGM\ 

N\m^\+SkT -^J+r'^qf-^r^Pr = F 


Here N - nu - constant, G is the gravitational 
constant, the mass of the sun, the radial 
component of the heat flux, and py. the radial 
component of the Poynting vector 

p = ^BX(uXB) 

In the right side of (5), F is a constant. Thus equation 
(5) expresses the conservation of the total energy in the 
steady-state expansion of the solar wind. During the 
expansion process, energy can be converted from one 
form into another, but the total energy must be 
conserved. 

The radial component of the Poynting vector is 
UjB^ sin^ <b 

Pr = - - 47 T 

Making use of (4), we can express the flow of magnetic 
field energy per steradian as 

r^Py. oc ^ (6) 

The solar wind velocity u is an increasing function of the 
heliocentric distance r. Thus the flow of magnetic field 
energy per steradian is a decreasing function of r. 
Conservation of the total energy flow requires that the 
field energy be transformed into other forms of energy 
flow. 

Mechanism of Energy Conversion 

The equation of motion for the magnetohydrodynamic 
expansion of the solar wind is 

“ m^-ng - V {InkT) + ^ JXB 

(7) 

where g = -eyGMJP . In the right hand side of (7), the 
magnetic field exerts a force 

i JXB 

on the moving solar wind flow. The radial component of 
this force can be written as 

sin^0 ^ 

4im dr 


( 5 ) This term is always positive, because u is an increasing 
function of r. The magnetic field exerts a force on the 
solar wind in the direction of the solar wind velocity. 


237 



The rate of work done by this force on a unit volume of ^ 
the solar wind flow is 


Mj 

t = -J-— 
^ 2kT* 


sin^ d) ^ 

47t dr 

The work done on the solar wind flow is used to increase 
the kinetic energy of the solar wind at the expense of 
the field energy . 


Because the denominator of (8) vanishes at the critical 
radius, we obtain 


1=14- sin^ 0* 


B*^ISn ^ 1 +2 sin^ 0* 
n^kT* 


In terms of the dimensionless quantities, we can write 
MAGNETOHYDRODYNAMIC ONE-FLUID MODEL the two governing equations, (8) and (9), in dimension- 
OF THE SOLAR WIND less form, 


Governing Equations 

The magnetohydrodynamical expansion of the one-fluid 
solar wind is governed by a system of two ordinary 
nonlinear differential equations : The equation of motion 

4kT _ _7MdT 

du^u 

- 2M. _ §1:^^ 

mi n 


dV 

dZ 


z 

z 


2d-'^-Z 


d§_ 

dZ 




1 


( 10 ) 


and 


M g cos^0* 
diZ~ cos^0 


Yl 

2 


+ (M)^ + 


i 

V 




( 11 ) 


and the equation of energy conservation 


K cos^ 0 


dT 

dr 


= N 


{rn,i 


4-nT’ 


r Aim 


)- 


F (9) 


where the thermal conductivity a is proportional to 

Let us denote conditions at the critical radius by a 
superscript asterisk. Then we can introduce the follow- 
ing dimensionless variables 


V^uju* 
6 = T/T^ 
Z = r/r* 


These two ordinary differential equations will suffice to 
solve for the two dependent variables F and 0 . 

At the critical radius, equation (11) reduces to 

= (12) 

Both the denominator and the numerator of equation 
(10) must vanish simultaneously at Z=l. When the 
numerator is zero, we obtain 

(i)'-2-2« (13) 

From the above two equations, we can show that a is 
not an independent parameter: 

“ " 4- K- 312) - 312 (14) 


and dimensionless parameters 

^ ^ - 2Nk . 
r*/c*cos^0'*' 

GMi^ 

H - F 
N m I 


Therefore, the general solutions for V and d can be 
expressed in the form 


and 


F= F(Z,//,7,|,0*) I 


(15) 


d = e(Z,H,y,l<PV 


Making use of the I’Hospital’s rule, we can calculate 
(dV/dZJ*-, the result is 
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where 


(iz = - Cl + (Cl " + C^f'^ (16) 


® + (l7 “ 2)(1 + 2 sin^0*) 
( 6 |- 2 ) 




2 + a(||'y-5^ + (^7 “ 2)^||7 + 2 sin^0* - 

(3|-1)~~ “ 


Numerical Solutions 

For any given value of % H, f , and 0 ^, we can use (13) 
and (16) to calculate the solution for V(Z) and in 
the e neighborhood of the singular point at Z = 1 . Then 
we can integrate equations (10) and (11) for Z > 1 + e 
and Z < 1 — 6 . The two branches of the numerical 
solutions join smoothly at the singular point. 

The numerical solutions for Z > 1 show the same 
property as the numerical solutions for the inviscid 
model of the solar wind [Whang and Chang, 1965] . As 
Z o®, the solution curve is expected to approach the 
condition F = Foo and 0= 0 with 



2 Woo 


A solution satisfying this condition when 7/ =3.0, 
% - 1.3, 0* = 176°, and 7 = 1.8978 is obtained as shown 
in figure 1 . The numerical values of this solution in 
dimensionless form F = V(Z) and 0 = d(Z) are given in 
table 1 . 

Physical interpretation of the numerical solution 
depends on the choice of the ratio rV''©* When the 
critical radius is chosen at between 5 and 6 , the solar 

wind conditions at 1 AU predicted by the present 
solution are (fig. 2) u = 287-3 18 km/sec and 
T = (1 .45 - 1 .57)X 10^ °K. In the remaining part of this 
paper, we will interpret the physical meaning of the 
solution based on r*/rQ=^5,5 (fig. 3). The predicted 
solar wind conditions at 1 AU are 
u = 302 km/sec 

r= 1.50X10^ °K 

^ = 5.7X10”® °K/cm 

(p = 129.5° 

B = 1.58 

where |8 is defined as the ratio of 2 nA:r to /Stt. 



Figure 1. The numerical solution, given in table 1, 
smoothly passes through the singular point at the critical 
radius (Z =l), and approaches the condition F=Foo 
an^ B = 0 as Z ^ 



Figure 2- The solar wind conditions at 1 AU predicted 
by the present solution depend on the choice of the 
ratio r^/r 0 . For r*jr^ = 5.5, the present result gives 
w= 302 kmisec and T= 1.50X10^ ""K at 1 AU. 
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Table 1 . Numerical solution for H- 3.0, t= 1 - 0 * ~ y = 1. 89 78 


z 

0 

V 

d 

deldZ 

& 

0.100 

171.3 

0.046 

3.026 

- 14.76 

0.021 

0.126 

172.7 

0.069 

2.704 

- 10.53 

0.020 

0.158 

174.0 

0.105 

2.415 

“ 7.483 

0.019 

0.200 

175.0 

0.159 

2.157 

- 5.294 

0.017 

0.251 

175.7 

0.234 

1.928 

- 3.733 

0.017 

0.316 

176.2 

0.331 

1.724 

- 2.630 

0.017 

0.398 

176.4 

0.448 

1.544 

- 1.853 

0.018 

0.501 

176.5 

0.579 

1.384 

- 1.309 

0.019 

0.631 

176.5 

0.719 

1.241 

- 0.926 

0.022 

* 0.794 

176.3 

0.860 

1.114 

- 0.657 

0.027 

1.000 

176.0 


1.000 

- 0.467 

0.033 

1.259 

175.6 

1.134 

0.898 

^ 0.333 

0.041 

1.585 

175.0 

1.261 

0.807 

- 0.237 

0.052 

1.995 

174.2 

1.380 

0.725 

- 0.169 

0.067 

2.512 

173.3 

1.489 

0.651 

- 0.121 

0.089 

3.162 

172.1 

1.590 

0.585 

- 8 . 61X10-2 

0.118 

3.981 

170.6 

1.682 

0.525 

- 6.14 1 

0.157 . 

5.012 

168.8 

1.766 

0.472 

- 4.38 

0.211 * 

6.310 

166.5 

1.842 

0.424 

- 3.11 

0.283 

7.944 

163.8 

1.911 

0.381 

- 2.22 

0.378 

10.00 

160.5 

1.974 

0.343 

- 1.58 

0.503 

12.59 

156.6 

2.030 

0.308 

- 1.13 

0.661 

15.85 

152.0 

2.081 

0.277 

- 8 . 10X10-2 

0.850 

19.95 

146.7 

2.127 

0.249 

- 5.88 

1.063 

25.12 

141.0 

2.169 

0.223 

- 4.35 

1.277 

31.62 

134.9 

2.207 

0.199 

- 3.27 

1.464 

39.81 

128.8 

2.241 

0.175 

- 2.51 

1.588 

50.12 

122.9 

2.271 

0.152 

- 1.96 

1.626 

63.10 

117.5 

2.298 

0.130 

- 1.53 

1.567 

79.44 

112.7 

2.322 

0.108 

- 1.18 

1.424 

100.0 

108.5 

2.342 

0.087 

- 8 . 78X10 -'* 

1.222 

125.9 

105.0 

2.359 

0.068 

- 6.16 

0.997 

158.5 

102.1 

2.372 

0.052 

- 3.87 

0.781 


The values of u, T, dT/dr, 0, and /3 calculated above are 
independent of the constant K used in the expression for 
the thermal conductivity k ~ . However, the 

proton number density, the magnitude of the mag- 
netic field and the heat flux at 1 AU predicted by 


the present solution depend on the value of the con- 
stant K (table 2). When the value of K is about half of 
that used in the ordinary expression for the thermal 
conductivity of fully ionized hydrogen, the predicted 
n, B, and q agree very well with the observed solar 
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TEMPERTURE CK) 




Figure 3. The solar wind temperature, velocity, direc- 
tion angle, and Rvalue as functions of the heliocentric 
distance, at I AU u = 302 km/sec, T= 1.50X10^ ^K, 
<t> = 129y and 


Table 2, The proton density, the magnitude of the 
magnetic field, and the heat flux at 1 AU predicted by 
the present solution 


ergs cm“^ 
sec“^ deg"'^^^ 

n, 

protons 

cm~3 

B. 

gammas 

ergs cm-* sec-' 

6.00X10-'^ 

13.1 

9.3 

2.23X1 0-* 

4.58X10-'' 

10.0 

8.1 

1. 70X10-* 

4.00X10-'' 

8.7 

7.6 

1. 49X10-* 

3.66X10-'' 

8.0 

13 

1. 36X10-* 

2.75X10-'' 

6.0 

6.3 

1. 02X10-* 


wind conditions [Hundhausen, 1970; Ness, 1967]. 
For K == 3.66X10“^ erg cm"^ see“^ deg“\ the 
predicted results at 1 AU are « = 8 protons/cm^ , 
B-13j and ^ = 136X10“^ erg/cm^ /sec. The pre- 
dicted n, B, and q as functions of the heliocentric 
distance r are plotted in figures 4 and 5. The density 



Figure 4. The proton density and the magnitude of 
the magnetic field as functions of the heliocentric 
distance. 



Figure 5, The heat flux as a function of the heliocen- 
tric distance, q = lAX lOr'^ ergjcm^ /sec at 1 AU. 
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curve agrees reasonably well with the observed data in 
the region 2 <r <20 [van de Hulst, 1950, 1953; 
Michard, 1954; Blacicwell, 1956] . In the region r< 2 
the present solution is not expected to be accurate, 
because the structure of the solar magnetic field is more 
complicated than the simple model described by 
equations (3) and (4). 

CONCLUDING REMARKS 

The solar wind conditions predicted by the present 
magnetohydrodynamic model show that when the effect 
of the solar magnetic field is included in the one-fluid 
model, we obtain an increase of 17 percent in the solar 
wind velocity or 37 percent in the convective kinetic 
energy at 1 AU. The variations of the magnetic field 
energy flow and the convective kinetic energy flow per 
steradian as functions of the heliocentric distance r are 
plotted in figure 6. The magnetic field energy flow per 



Figure 6. The radial flow of magnetic field energy per 
steradian is a decreasing function of the heliocentric 
distance. The field energy is converted into the kinetic 
energy during the expansion of the solar wind. The 
convective kinetic energy flow per steradian at 1 AU is 
chosen as unity in this figure. 


steradian is a monotonically decreasing function of r. 
Using the convective kinetic energy per steradian at 
1 AU as a unit, we can see that the total decrease of the 
magnetic field energy flow from r-lr^to r=\ AU is 
about 0.35 . Conversion of the magnetic field energy into 
the kinetic energy certainly plays an important role in 
the expansion process of the solar wind. 

Urch [1969] has studied the magnetohydrodynamic 
one-fluid model of the solar wind. The equations of 
motion and energy conservation were integrated in a 
quite different way. He has obtained the numerical 
solutions for the solar wind after several iterations. The 
temperatures at 1 AU predicted by Urch are 3 to 
4X10^°K, which are too high compared with the 
observed data. 

Recently Wolff et al [1971] have studied the 
magnetohydrodynamic two-fluid model of the solar 
wind; their predicted solar wind conditions at 1 AU are 
very close to the results obtained in this paper. Barnes et 
al [ 1971] have also included the magnetic field in their 
two-fluid model of the solar wind. Because they have 
dropped out the JXB/c force tenn in their equation of 
motion, the magnetic field force is not doing work on 
the expanding solar wind. The conversion of magnetic 
field energy into kinetic energy is not included in their 
work. They proposed that an energy source due to 
dissipation of hydromagnetic waves supplies the energy 
needed to increase the solar wind velocity. Their source 
of energy is completely different from the magnetic field 
energy flow (the Poyn ting vector) discussed here. 
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